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NON PERTURBATIVE EFFECTS IN QCD 

Adriano Di Giacomo^ 

"'Dipartimento di Fisica Universita' Pisa and INFN Sezione di Pisa, PISA ITALY 

Non perturbative results from lattice QCD will be discussed, namely : Vacuum Condensates and QCD Sum Rules; (7a (1) 
and Topology; Confinement of Color. 



1 INTRODUCTION 

QCD is usually quantized perturbatively. 

The Lagrangean L is split into the sum of two terms 

L = Lo + Li (1) 

where Lq describes free quarks and gluons ,and Lj 
their interaction. Then Lq is quantized : the Hilbert 
space is the Fock space of free quarks and gluons and 
the interaction describes scattering between them. 

Any observable (O) is computed as a power series ex- 
pansion in the renormalized coupling constant as 

(O) = S„0„a^ (2) 

where On are finite amplitudes. 

The expansion eq(2) is not convergent, not even as an 
asymptotic series[^. Nevertheless it is an empirical 
fact that, as — > oo or for > IGev^ as ^ 1 and a 
few terms of the expansion describe physics correctly. 
The lack of convergence of the expansion eq(2) reflects 
the instability of Fock vacuum. Quarks and gluons are 
confined and therefore free quarks and gluons are not 
a good zeroth order approximation to physics. 

A non perturbative quantization of the theory is 
needed to identify the true vacuum state. This is pro- 
vided by the Feynman path integral formulation. The 
key quantity is the partition function Z, in terms of 
which all field correlators can be computed. 

Z = J[dA^] [dV] m exp - [S{A^,{p, V-)] (3) 

Z in eq(4) is a functional integral and is defined by 
approximating euclidean space time by a discrete set 
of points in a finite volume, so that the number of 
integration variables is finite and the integral is an or- 
dinary well defined integral. A sequence of such inte- 
grals is then computed on sets of points which tend to 
cover the volume densely. The limit of infinite volume 



is then performed. If Z is finite and well defined after 
these limiting procedures one says that the theory ex- 
sists as a field theory. The procedure also ensures the 
existence of the analytic continuation to Minkowskian 
space time. 

For QCD asymptotic freedom insures that the first of 
these limits exists. The existence of a mass gap in 
the theory also insures that the infinite volume limit 
exists . A fully rigorous mathematical proof of these 
statements still does not exist, but it is physically rea- 
sonable to say that QCD exists as a field theory. 

Lattice formulation of QCD is nothing but an approx- 
imant of Feynman integral in the sequence which de- 
fines it : at sufficiently low value of the unrenormalized 
coupling and at sufficiently large volume it will provide 
the physical amplitudes with any required precision 
from first principles (modulo technical computational 
difficulties) . 

Non perturbative phenomenology is based on Wilson's 
Operator Product Expansion (OPE) , 

TiA{x)B{0)) = I]„C„(a;)O„(0) (4) 

where C„'s describe short distances and are usually 
computed in perturbation theory and the operators 
On describe large distances. The sum in eq(4) is or- 
dered by incresing order of the dimensions in mass of 
the local operators 0„. 

Vacuum expectation values (vev) of eq(4) involve in 
the rhs vev's of composite operators (condensates) . 
They enter in the SVZ sum rules The gluon con- 
densate (G2) = (0|/3(g,)/.g,G^,,(0)G^,(0)|0) , and the 
quark condensate (G^^) = {0\ip{0)tp{0)\0) are exam- 
ples of them. A condensate like G2 has dimension 4 in 
mass, and therefore, by renormalization group argu- 
ments is related to the running coupling constant at 
the scale /i as G2 = /i^ea;p(— 4/6o/5^(a*^)) which is non 
analytic in g and undefined in perturbation theory. 

Matrix elements of eq(4) between quark states are rel- 
evant to weak interaction matrix elements at distances 
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= l/Mw^ 1/A/z ,or to structure functions of deep in- 
elastic scattering. 

In this talk I will concentrate on the vacuum state, 
namely on condensates (sect 2) , topology (sectS) , 
and confinement of color (sect4) 

2 QCD SUM RULES. CONDEN- 
SATES 

The OPE of the product of two conserved currents 

r(jM(a:)>(0) = (g^. - d^d,)Ii{x^) (5) 

taken on the vacuum reads 

{tl{x^)) ^ C}{x^m+CG{x^){G2)+C^^{x^){G^^)+..{&) 

where I is the identity operator ( (/) = l),and at 
short distances C/(a;^) ^ l/x** modulo logs,CG and 
~ const modulo logs. Eq(6) is a theorem in per- 
turbation theory, an assumption in the presence of 
non perturbative effects. The coefficients C describe 
short distances, the condensates describe large dis- 
tance physics. 

After Fourier transform 



n(g2) - n(0) = j d'^x{n{x'^)){exp{iqx) - 1) ~ 
Ci{q^) + CGAq'){G2) + C^4q'){G^^) (7) 

with 

Gi{q^) ^ const (modulo logs) (8) 
Cg.('z'), C^-.4q^) ^ 1/q^ (modulo logs) (9) 

A dispersive representation for the l.h.s. of eq(7), if 
is the electromagnetic current is 



n(g2)_n(0) = -g2y dfi' 



i?(/i2) 



//^(g^ - /i^ + ie) 



(10) 



where R{s) = ' " can be taken from exper- 

iment. 



An appropriate average on q^ of both members of 
eq(lO) gives the SVZ[|2] sum rules, in which non per- 
turbative effects are parametrized in terms of the con- 
densates {G2) and (G^^). The result is a good phe- 
nomenology and a determination of the condensates 



(G2> - (.024±.011)Ge'i;^ 
{G^J -.13Gew3(g2 = iGev^) 



(11) 
(12) 



However the perturbative expansion of Hiq"^) — 11(0) 
which provides the coefhcient Ci{q^), when resummed 
at higher orders , is ambiguous by terms oc 1 [1 , [ 
|S] which mimic the terms with G2 and G^^. The defi- 
nition of the condensates is then intrinsically ambigu- 
ous. They could be defined as the coefficients of the 
terms 1/q^ in the expansion eq(7), but a priori, due to 
the ambiguity, they could be process dependent, i.e. 
dependent on the currents under consideration. 

An OPE of the levels of bound QQ energy levels [0] [IH] 
can be also performed. This amounts to compute the 
quadratic Stark effect produced on the levels by the 
fluctuating vacuum chromoelcctric field. If the finite 
correlation length of the field strength correlators is 
taken into account [ E] [ d the relevant quantity is the 
gauge invariant correlator [|H] 

{G2{x)) = {0\Tr{E{x)Uc{x,0)E{0)U\x,0)}\0) (13) 

with Uc{x, 0) the parallel transport from to a; along 
the path C, or Uc{x,Q) = Pexpii ^ ^ A^{x)dx'^) . A 
general parametrization of such correlators is [El [EI] 

D^,,pa{x) = {Q\Tr[G^,,{x)Uc{x,Q)Gp„{{))ul:{x,Q)m - 
{9t.p9i^<y - 9t.<y9up)[D{x^) + Di{x^)] + (14) 

~t~ fj,'^ (7 9 !^ p •^L''^a9fJ^p '^U'^ p9 p<t\ Qx^ 



In this language (G2(x)) = aj-nip -f D\). The corre- 
lators Dpypa can be computed on the lattice [I5][[TT] 
and with them their OPE. Their typical behaviour is 
shown in figs 1 and 2. A convenient parametrization 
of the Lattice data D^, is [Ellini 

-^D"^ ^^ + C7G(G2)exp(-a:/A) (15) 
a* x^Q X* 



^^^ + G'G(G2)exp(-x/A') (16) 



{G^ix)) is a split point regulator of (G2). A 
best square fit to the lattice data gives, for pure 
gauge SU(3)[|E| (quenched, no quarks) G2 = (.15 ± 
.03)Ge^;^ A = A' = (.22 ± .03)/m . 

In full QCD [ ^] the same quantities can be com- 
puted. An extrapolation to realistic quark masses 
gives (G2) = (.022 ± .005)Gew'* and A = A' = 
(.32 ± .04)/m. The phenomenological value for (G2) 
is (.024± .005Gew)4[[T3] 

From quenched to full QCD the gluon condensate de- 
creases by almost an order of magnitude, and the cor- 
relation length increases by 50%. 

Lattice vacuum correlators are the input of the 
Stochastic approach to QCD [CHI [EI- 
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Figure 1. D^{x). The line is the best fit to eq(15). 
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Figure 2. D^{x)i. The line is the best fit to eq(16). 



In a similar way the fcrmion correlators can be com- 
Duted[rT5| 

{G^.4x)) ^ (0|^(x)[/c(:c, 0)^(0)10) (17) 
{G^jp{x)) can be viewed as a split point regulator of 

Lattice data [ E| give {0\ipilj\0) consistent with the 
value obtained by use of SVZ sum rules and a corre- 
lation length A = (.42 ± .05) /m 

3 TOPOLOGY 

The 1/^(1) problem was a problem of the free quark 
model of Gellmann. The singlet axial current is con- 
served in the model, but in Nature the corresponding 
symmetry is neither realized a la Wigner (no parity 
doublets exist in the hadron spectrum), nor a la Gold- 
stone : indeed a Goldstone-broken symmetry would 
imply [E| w,,' < \/3™7r , which is badly violated by 
the observed value m^/ ^ 980Meti. Hence the idea 
of Gellmann that the symmetry of hadrons could be 
abstracted from the free quark model was not correct 
for the axial singlet. 

QCD solves the puzzle. The singlet axial current is 
anomalous in QCD and UAi^) is not a symmetry,being 
broken at the quantum level. 

d,,r{x) - 2NfQ{x) (18) 

(G*j, = l/'2.e^yp„Gpcr ) is named topological charge 
density. 

On continuous configurations the topological charge 
Q — j d^xQ{x) is an integer (second Chern Number). 

For quantized fields Q{x) is an operator. The Topo- 
logical susceptibility x is defined as the response of the 
vacuum to Q(x) 

X = Jd-'xK 0|r(Q(x)Q(0))|0 > (19) 

An argument based on Nc — s- oo shows that m^/ is 
related to the topological susceptibility of the Nc — oo 
vacuum. [IT^ 

2NfX = fW,'a + 0a/Nc)) (20) 

As Nc — !■ oo dynamical quarks can be neglected and 
with the same approximation Nc can be taken equal 
to the physical value 3 ,so that x is the topological 
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Figure 3. Determinations of x by use of different lattice 
regulators and of different values of the unrenormalized 
coupling constant. The physical value of x is independent 
of these choices within statistical errors. 



susceptibility of the quenched SU(3) vacuum. An im- 
proved version of the argument including mixing to 
the octet gives [H^ 



2NfX = fUK 2ml)il + 0(l/iV,)) (21) 

The prediction is in both cases x ^ (ISOMew)^. 

Computing x from first principles can provide a cross 
check of QCD and of the l/Nc expansion at the same 
time. 

This problem has a long story: tha main difficulty 
has been for a long time to understand how topology, 
which is based on continuity, can exist on the lattice 
which has a discrete structure. The way out is basic 
field theorv[On| : any lattice version of the operator 
Q{x), Ql{x), is related to the continuum counterpart 
,say in the MS scheme, by a multiplicative renormal- 
ization 



Xx) = ZqQ{x) 



(22) 



The product Q(x)Q(0) is singular as a; O.The OPE 
gives 

Q(x)Q(0) - ci{x)I + C2{x)G2 + finite terms (23) 

„-4 



Cl ^ X ,C2^X 



(24) 



to compute Zq and Xl is known as " heating" [ : it 
is based on the fact that topology has a much longer 
autocorrelation time tq in montecarlo upgrading than 
the autocorrelation time r of the local quantum fluctu- 
ations which dominate the renormalization constants. 
Therefore can be measured on the configurations 
obtained at times t such that t <ti t ^ tq from a 
zero field configuration, which has topological charge 
zero. Similarily Ql can be measured on configura- 
tions prepared from an initial configuration consisting 
of a single instanton , which has Q=l. At times long 
enough with respect to r but very short compared to 
Tq Ql will be equal to Z . This method is known 
as "field-theoretical". The result of this procedure is 
X — ((175 ± 5)Mev)'^[ 24^ Recently an independent 
determination has been made, using the divergence of 
the singlet axial current operator, i.e. the Ihs of eq 
(17). The advantage is in the use of a formulation 
of the fermions on the lattice which preserves chiral 
symmetry, and hence has multiplicative renormaliza- 
tion Zq = 1 The result is[[2K| 



X = ((190± 10)^/6^) 



(26) 



In computing x additive renormalization is also 
required. [ 21 In conclusion [I22| 



in complete agreement with the field theoretical deter- 
mination. Also determinations based on cooling tech- 
niques [ 121] give consistent estimates. In conclusion 
the C/^(l) problem is solved in QCD, and additional 
legitimation is added to the ideas of A^c ^ oo. An 
important issue, as we shall see in the next section, is 
the behaviour of x at the deconfining phase transition. 



4 CONFINEMENT OF COLOR 

Quarks and gluons are visible at short distances, but 
have never been observed as free particles. This fact 
has lead to the conjecture that colored particles never 
appear in asymptotic states, a property known as Con- 
finement of Color. Confinement should be derivable 
from the QCD lagrangean. 

Quarks have been searched in Nature since they were 
introduced by Gellman as fundamental constituents of 
hadrons [ i27j , with negative results: only upper limits 
to their abundance and to their production rate have 
been established [1^. 

The ratio of quark abundance in nature, Uq to that of 
protons Up has been given the upper limit 



X = 



X L - Xl 

Z2 



(25) 



where Xl is the lattice topological susceptibility of the 
Q ~ Q sector of the vacuum. A technique developed 



< 10 



-27 



(27) 



to be compared to the expectation based on the Stan- 
dard Cosmological Model ^ w 10" . 
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As for production in particle reactions the best limit is 
provided by the inclusive cross section Uq = (j(j)+p ^ 
q{q) + X) namely 



ag < 10-42.^2 



cm 



(28) 



to be compared with the expectation, in the absence 
of confinement, 



CTo — CTtOTAL ~ 10 ^^CTO^ 



(29) 



In both cases the ratio is depressed by a factor ~ 
10^^^. The only natural explanation is that these ra- 
tios are exactly zero, i.e. that confinement is an abso- 
lute property, and as such based on some symmetry. 

Experiments have been and are being performed to 
detect a deconfining transition at high energy density 
from hadronic matter to a state of free quarks and glu- 
ons (Quark Gluon Plasma), in high energy collisions 
of heavy ions [ 123 • The major difficulty with them is 
to identify an observable quantity by which deconfinc- 
ment could be detected. 

Some evidence for the existence of that transition has 
been produced by virtual experiments, namely by nu- 
merical simulations of QCD on a lattice. 

The partition function of a field theory at temperature 
T is the euclidean Feynman path integral extending 
in euclidean time from r = to r = y, with periodic 
boundary conditions for boson fields, antiperiodic for 
fcrmions. 

For QCD this is computed on a lattice of spatial ex- 
tension and temporal extension Nt with Nt Ng, 
and the temperature is T = ^^^^^^^^ with a(/?)the lattice 
spacing in physical units , which can be tuned by vary- 
ing the unrenormalized coupling constant /? 



Also on the lattice, however, as in experiments, the 
real problem is to have a criterion to detect confine- 
ment. 

In quenched formulation (pure gauge, no dynamical 
quarks) the criterion consists in looking at the large 
distance behaviour of the static QQ potential, V{x) 
. V(x) is related to the correlator of Polyakov lines 
Dix) 



D(x) 



(L(f)it(0)) 



by the relation 

V{x) = ~^bi{D{x)) 



(30) 



(31) 



where aNf is the extension of the lattice in the time 
direction. By definition the Polyakov line is the trace 



of the parallel transport across the lattice in the time 
direction. 



L{x)=Pexp[i dTA°{x,T)] 
Jo 

By cluster property 

D{x) ~ exp{-ax/T) + \{L)\^ 

X— i-OO 

If (L) = it follows 



(32) 



(33) 



V{x) ^ ax as a; ^ oo (confinement) (34) 
If {L) ^ 

V{x) ~ const, as x oo (deconfinement) (35) 

A transition is observed on the lattice at Tc ~ 270M ev 
from a region T < Tc where {L)=0 to a region T > Tc 
where (L) ^ 0. A finite size scaling analysis allows to 
determine the order of the transition, which for SU(3) 
is first [ 1^ [ 1201 • (L) can be then assumed as an or- 
der parameter for confinement , and the corresponding 
symmetry is Z3 

In principle one should show that (L) — implies 
absence of colored particles in all asymptotic states to 
prove confinement, which has not been done, but the 
criterion is reasonable anyhow. 

In the presence of dynamical quarks (full QCD) Z3 
is not a symmetry anymore and (L) is not an order 
parameter. String breaking is expected to take place: 
when pulling the static QQ pair apart from each other 
energy is converted into light qq pairs, and the poten- 
tial is not linear any more, even if there is confinement. 

At niq — another symmetry is present, chiral sym- 
metry, which is spontaneously broken at T = , and, 
as lattice data show, is restored at some Tc- Chi- 
ral symmetry, however, is explicitely broken by quark 
masses, and cannot be the symmetry responsible for 
confinement. 

Let us consider the case Nf = 2, m„ = nid ~ m, which 
is semirealistic, but, as we shall see, very instructive. 

A phase dyagram is usually drawn for this system as 
in fig(4). The line which joins the first order phase 
transition at m = 00 (quenched) to the chiral tran- 
sition at m = corresponds to the values Tc{m) 
at which the susceptibilities Cy (specific heat) and 
= J (Px{ip{x)i/j{x)il;{0)i/j{0)) are maximum [ '32' [ 



Conventionally the region below this line is called con- 
fined, the region above it deconfined. No criterion ex- 
ists which justifies this attribution. 
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Figure 4. Schematic phase diagram of Nf = 2 QCD. The 
line corresponds to maxima of the susceptibilities. 

An analysis can be made of the chiral phase transition 
at m = 0, assuming that the pions and sigmas are 
the relevant degrees of freedom[|^. Chiral symmetry 
fixes then the form of the free energy and renormaliza- 
tion group arguments allow to predict the order of the 
transition and the corresponding universality class. 

For Nf — 3 the transition is first order. For Np — 2 
the order depends on the relative weight of the chi- 
ral 0(4) part of the action and the term describing 
the C/4(l) anomaly. If the anomaly disappears below 
Tc the transition is first order, and such is the transi- 
tion along the line at m ^ 0. If instead the anomaly 
persists up to Tc the transition is second order the 
universality class is 0(4) and the line at m 7^ is a 
crossover. A reliable numerical analysis of this issue 
is not yet available, but for some reason the second 
option is more popular. 

The order and the universality class can be investi- 
gated by studying the behaviour of the critical spe- 
cific heat as the volume goes to infinity by a finite size 
scaling analysis. A good order parameter should have 
a behaviour in that limit consistent with that of Cy 
.The order of the transition and the universality class 
should be identified by the symmetry responsible for 
confinement. 

A valid candidate symmetry is dual superconductivity 
of the vacuum[ 35 [ 36i[|2Z|- The basic idea is that the 
chromoelectric field acting between colored particles is 
channeled into dual Abrikosov flux tubes, whose en- 
ergy is proportional to the length. Here dual means 
interchange of electric with magnetic with respect to 
ordinary superconductors. 

This mechanism has been assumed as a working 
hypothesis and analyzed down to observable conse- 
quences in a series of papers [ EH] [ EHI [ EOI [ EH [ El [ 

III][ESI- 

A disorder parameter (/i) has been defined which de- 
tects dual superconductivity: it is the vev of an op- 



erator fj, carrying magnetic charge which has zero vev 
in a phase in which vacuum has a definite magnetic 
charge, and can have a non zero vev if monopoles con- 
dense i.e. the vacuum is a superposition of states with 
different magnetic charge. 

That operator is color gauge invariant , magnetically 
charged but magnetic U(l) gauge invariant[ 11 . 

In principle (/i) ^ or(/x) — i.e. superconductivity 
or not depends on the procedure used to identify mag- 
netic charges (Abelian Projection). However it can be 
shown that monopolc condensation (or non condensa- 
tion) is an abelian projection independent statement [ 

E11[ES1- 

For quenched theory the result of numericall simula- 
tios is that this disorder parameter is consistent with 
(L)[EI1[E2 : it is non zero for T < and strictly 
zero for T > Tc- By finite size scaling analysis of the 
corresponding susceptibility p — j^log{{fi)) the same 
critical indices are obtained as with (L). 

The operator fi is perfectly well defined also in the 
presence of dynamical quarks, and can be used to in- 
vestigate the phase diagram of fig(4). The result is 
that the region below the line is really confined, the 
one above it is deconfined[l46j. 

p as a function of T at fixed m has a peak at the 
same value as the other susceptibilities, i.e. on the 
line of fig(4). (/i) thus provides a good criterion for 
confinement. 

The finite size scaling analysis works as follows. A 
priori (/i) depends on the coupling constant /? i.e. on 
the lattice spacing a , on the mass m and on the lattice 
size Ns .The dependence on /? can be traded with the 
dependence on the correlation length ^ so that, for 
dimensional reasons 

(/i) = $(|,^,^7V7) (36) 

As T approaches Tc from below, or t = (1 — ^) — > 0, 
diverges as 

^ cx T'" (37) 
v and 7 are critical indices . 

If the transition is second order or weak first order ^ 
goes large as Tc is approached and | can be put equal 

to zero. Because of eq(37) ^ can be traded with 
tNs^'^ ■ The scaling law follows 

p/7Vy^ = /(rA^y^m7V^) (38) 

In the quenched case the second dependence does not 
exist. The scaling can be tested on the lattice data and 
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Figure 5. The behaviour of the topological susceptibil- 
ity (cercles), of the magnetic monopole susceptibility p 
(squares), and of the chiral parameter {iji'ip) (triangles) ver- 
sus r/Tc at the deconfining transition. 

the critical index v can be extracted[E][02|- In Par- 
ticular for the peak one finds 

Ppeak (X N]/'' (39) 

For the case of full QCD the problem has two scales 
and the analysis is more complicated. Preliminary re- 
sults[23 exclude the critical indices of 0(4) and seem 
consistent with a first order phase transition. If this 
will be confirmed by the additional data in prepara- 
tion either the anomaly disappears below , or the 
relevant degrees of freedom which dominate are not 
the pion and the sigma[IS^. 

The first possibility can be checked on the lattice. The 
dependence of the topological susceptibility at the de- 
confining transition has been studied[|^.Fig(5) shows 
the existing data. A more precise analysis is needed 
to settle the problem. 

In conclusion there is strong evidence that dual su- 
perconductivity is the mechanism of confinement. A 
candidate disorder parameter exists , which is well 
defined and seems to work. The corresponding sym- 
metry is dual superconductivity, which is however not 
clearly understood formally, being a common feature 
of all the abelian projections: consequently the effec- 
tive action is not known. 

A definite settlement of the order of the transition for 
the case TV/ — 2 will be a crucial test. 
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